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Abstract 

The Dirac equation in a curved spacetime depends on a field of co- 
efficients (essentially the Dirac matrices), for which a continuum of 
different choices are possible. We study the conditions under which 
a change of the coefficient fields leads to an equivalent Hamiltonian 
operator H, or to an equivalent energy operator E. We do that for 
the standard version of the gravitational Dirac equation, and for two 
alternative equations based on the tensor representation of the Dirac 
fields. The latter equations may be defined when the spacetime is 
four-dimensional, noncompact, and admits a spinor structure. We 
find that, for each among the three versions of the equation, the vast 
majority of the possible coefficient changes do not lead to an equiv- 
alent operator H, nor to an equivalent operator E, whence a lack of 
uniqueness. In particular, we prove that the Dirac energy spectrum 
is not unique. This non-uniqueness of the energy spectrum comes 
from an effect of the choice of coefficients, and applies in any given 
coordinates. 

1 Introduction 

Dirac's equation for a relativistic spin-half particle was originally proposed 
for a Minkowski spacetime. Dirac's original equation is limited to Carte- 
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sian coordinates, and the coefficients which enter that equation — the Dirac 
matrices — are constant. In a curved spacetime, or for non-Cartesian coor- 
dinates in a Minkowski spacetime, Dirac's original equation does not apply 
any more, and is usually replaced [H 121 [3] by the equation proposed indepen- 
dently by Weyl [1] and by Fock [5] , hereafter the "Dirac- Fock-Weyl" equation 
[DFW for short). This equation has been investigated in situations of physi- 
cal interest, e.g. in uniformly rotating coordinates in a Minkowski spacetime 
(e.f?. [3]), in uniformly accelerating coordinates in a Minkowski spacetime or 
in a static, or stationary, weak gravitational field {e.g. [21 [6l [TJ El HI [TO]), 
in a Robertson- Walker spacetime [TTl [T2] . etc. While the currently observ- 
able effects of the gravitational field on quantum particles (such as neutrons 
or atoms) can still be described by using the non-relativistic Schrodinger 
equation in the Newtonian gravity potential, one may expect that, in a not- 
too-distant future, improvements in the experimental accuracy should allow 
one to check the combined effects of relativistic gravity and relativistic quan- 
tum mechanics. This will be a check of the way in which we conceive the 
coupling between gravity and the quantum, including a check of the way we 
write the quantum wave equations in a curved spacetime. 

The foregoing references [1]-[I2] address the first quantized Dirac theory. 
Indeed, there is a long list of physical effects whose interpretation uses the 
first quantized Dirac theory and does not require second quantization — e.g., 
the Stern-Gerlach, COW, Sagnac, and Mashhoon effects. Moreover, it is 
well known that the problems of quantum field theory in a curved spacetime 
do not all reside within quantum field theory itself, since some of the most 
challenging problems for quantum fields in a curved spacetime involve diffi- 
culties in defining the single particle states, along with their Hilbert space 
structure, and the operators that act on them. These difficulties with sin- 
gle particle states, which are used to define the field operators, the vacuum, 
and the occupation states, are not resolved by the mere formalism of quan- 
tum field theory. Among these are the difficulties associated with the energy 
spectrum of single particle states in a curved spacetime. This is critical, not 
only for its impact on proposed experiments, but also theoretically for the 
identification of positive and negative energy states. See for example Refs. 
[131 HI]- Therefore, just like Refs. [1]-[12], the present paper addresses the 
first quantized Dirac theory, which gives the normalized single particle and 
antiparticle states that define a second quantized fermion field operator in a 
curved spacetime [121 [E] . 
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In the DFW equation, the relevant, "curved" Dirac matrices 'j'^{X) most 
generally vary from one spacetime point X to another one, even though they 
are expressed linearly as functions of a fixed set of constant Dirac matrices, 
say 7''", through a matrix of coefficients a^^{X) associated with an orthonor- 
mal tetrad field: 

nX) = a^^iX) 7^". (1) 

In a given Lorentzian spacetime {Sig^ip)-, there is a continuum of possible 
choices for the tetrad field and its matrix a{X) = {a'^a) — ^^is has been em- 
phasized recently by other authors also |15l |16]. The tetrad field is assigned 
to satisfy the orthonormality condition: 

a^,a'^^r/"'^ = ^7^^ {g^'^) ^ {g^,)-\ (2) 

[Here (r/"^) = {ria/3)~^ = {rjap) = diag(l, — 1, — 1, — 1).] Moreover, any set 
(7^") of "fiat" Dirac matrices [i.e., any constant solution of the anticommu- 
tation relation (|4]) below with g'^'^ = t]^'^] enjoys exactly the same status. 
Hence, none should be preferred over another one [17], although in the liter- 
ature the choice of the set (7''") is restricted unnecessarily, as we will see in 
Section \2.A\ by a hermiticity condition [181 [H] • 

As stated above, in Eq. ([T]) there is a continuum of possible choices for 
the field 7^. This is also true for the two alternative extensions of the Dirac 
equation to a curved spacetime, which have been proposed recently [20] . 
Each of these two alternative Dirac equations in a curved spacetime is based 
on the tensor representation of the Dirac (TRD) field [201 EI], in which the 
wave function is a four-vector and the set of the components of the four 7^^ 
matrices builds a third-order tensor. 



^ In a Minkowski spacetime in Cartesian (or affine) coordinates, the TRD theory with 
constant Dirac matrices has been proved [17 to be quantum- mechanically fully equivalent 
to the genuine Dirac theory. See also Subsect. 1.1 in Ref. [20] for a summary of the 
argument. Like the solutions of the genuine Dirac equation, the single particle TRD 
solutions have only two spin polarizations (up and down) which makes them spin-half 
wave functions. Moreover, just as with the genuine Dirac equation, it is straightforward 
to extend the single particle theory to a canonical second quantized fcrmion theory (see 
Subsect. 4.3 of Ref. [T^), with the fermion field operator built on the normalized single 
particle and antiparticle TRD solutions. 
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In a Minkowski spacetime, the Dirac wave function transforms under 
the spin group as the coordinates are Lorentz-transformed. However, it is 
well known since the work of Fock and Weyl [H |5] that, on changing the 
coordinates, transforming the Dirac wave function under the spin group is 
not an option in a curved spacetime or even in a flat spacetime with affine 
coordinates. Thus, there are only two possibilities in a curved spacetime. 
Either the Dirac wave function transforms as a quadruplet of four scalar 
fields under coordinate transformations as in DFW, or it transforms as a 
complex four- vector field as in TRD. As discussed in Subsect. 12. the lat- 
ter four-vector transformation is possible whenever the spacetime manifold 
is four- dimensional, noncompact, and admits a spinor structure. In a flat 
spacetime with affine coordinates, and constant Dirac matrices, TRD and 
DFW have been shown to be equivalent [VT}. From the Lagrangian in Sub- 
sect. 12. 3[ we can show in a straightforward manner that this equivalence 
gives rise to equivalent fermion quantum field theories in a fiat spacetime 
with affine coordinates. 

However, to our knowledge, there has been hardly any study of the pos- 
sible effect of the choice of the field 7^^, even for the standard, DFW theory. 
Moreover, when one uses a general set of Dirac matrices, one needs a matrix 
A, called the hermitizing matrix, that also becomes a field of matrices A{X) 
in a curved spacetime. So the Dirac equation in a given spacetime depends on 
the fields (7'^, A), which we call the coefficient fields of the Dirac equation. In 
a recent work [22] , the quantum mechanics associated with these three grav- 
itational Dirac equations — the standard, DFW equation, and the two TRD 
equations — has been studied with regard to conservation of the probability 
current, Hilbert space scalar product, and hermiticity of the Hamiltonian op- 
erator. In particular, the hermiticity condition has been derived in a general 
spacetime and it has been shown that, for DFW, the hermiticity of the Dirac 
Hamiltonian is not preserved under all admissible changes of the coefficient 
fields 

This surprising result indicates that the standard version of the gravita- 
tional Dirac equation has a non-uniqueness problem. The aim of the present 
work is, therefore, to begin the systematic study of the non-uniqueness prob- 
lem of the Dirac theory in a curved spacetime. The result just recalled points 
to a non-uniqueness of the Hamiltonian theory, thus, possibly, of the energy 
spectrum associated with the Dirac equation in a curved spacetime. In a 
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time-dependent metric, the Hamiltonian operator, that is the generator of 
the time evolution, is not necessarily Hermitian and does not generally co- 
incide with the energy operator, associated with the field energy. It is the 
energy operator, not the Hamiltonian, which is derived canonically from the 
Lagrangian: see Ref. [23], and Sect. H] below. Hence, in the present work, 
we focus on the non-uniqueness problem of the Hamiltonian and energy op- 
erators, both having physical relevance. We ask successively: When does a 
change of the coefficient fields lead to an equivalent Hamiltonian operator? 
Or to an equivalent energy operator? In any given coordinate system is the 
energy spectrum unique? We shall investigate these questions for the DFW 
theory as well as for the alternative, TRD, theory. In particular, in Sect. 
we prove that, in any given coordinate system, the spectrum of the en- 
ergy operator for the free Dirac equation is not unique in both DFW and 
TRD, but depends on the tetrad field or set of coefficient fields that is chosen. 

While it is expected that the energy spectrum would depend on the choice 
of the time coordinate, it is unexpected to discover that, in any given coor- 
dinate system (in particular, with any given time coordinate), the energy 
spectrum varies with the choice of the coefficient fields {'-f^,A). Currently, 
there is no preferred way to choose these coefficient fields in a general curved 
spacetime. The non-uniqueness of the energy spectrum is significant be- 
cause it is a fundamental axiom in quantum mechanics that the spectrum of 
the energy operator at a given time t restricts the allowed values of energy 
that will be observed in experiments at time t. This axiom, first enunciated 
by Max Born, is postulated to be true of both time- dependent and time- 
independent quantum theories. Consequences of Horn's axiom together with 
the non-uniqueness of the energy spectrum in any given coordinate system 
are discussed briefiy at the end of Sect. Ei We show that the energy shifts 
due to different "reasonable" choices of the coefficient fields (7^^, A) are very 
small, similar in magnitude to the Mashhoon energy term for a rotating coor- 
dinate system in a Minkowski spacetime. Such small energy shifts, described 
in Sect. El are too small to observe currently, but may be detected in the 
future with improved experimental techniques [21] 
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2 Hamiltonians of the Dirac equation and the 
modified Dirac equation 

2.1 Hamiltonian of the ("normal") Dirac equation 

As is well known, a given spacetime {y-,g^u) need not admit a spinor struc- 
ture. However, the situation is greatly simplified if the spacetime is four- 
dimensional and noncompact. It was proved by Geroch that such a spacetime 
admits a spinor structure if and only if it admits a global tetrad field |25] . 
In that case, both spinor bundles and tangent bundles are trivial [251 126] . 
Penrose and Rindler argue that these are the only spacetimes of interest [27] . 
The normal Dirac equations in a curved spacetime (V,(7^j,), either the stan- 
dard, DFW equation, or any of the two alternative TRD equations, have the 
same form: 

7^D^V^ = -imi), (3) 

where 7^ = 7^(X) (yU = 0, 3) is a field of 4 x 4 complex matrices defined 
on the spacetime manifold V, satisfying the anticommmutation relation 

7V + 7V = 2(/'^n4, /i,J/e{0,...,3} (l4 = diag(l, 1,1,1)); (4) 

and where m is the rest-mass of the particle (setting h = 1 = c). In a 
curved spacetime, if) is either a quadruplet of four scalar fields (for DFW) or 
a four- vector field (for TRD). Thus, for DFW, the field if) remains invariant 
after a coordinate change IH [21 13] and, accordingly, the 7^^ 's transform then 
as a four-vector. Instead, for TRD, the field ip is a. true four-vector, and 
the array of components 7^'' = {■j'^Y^, builds a (J) tensor [201 If there 
is a global tetrad field on the spacetime V (as proved in Geroch's theorem 
recalled above), then it is easy to show that one may define a global field of 
Dirac matrices 7^^ satisfying Eq. (jl]) above, for DFW and for TRD as well. 
The three equations differ also in the connection that defines the covariant 
derivative D^. For each of the two TRD equations, this is a connection on 
the complexification TcV of the tangent bundle — namely, the Levi-Civita 
connection extended to TcV for "TRD-1", and a special affine connection 
for "TRD-2" [201 [22]. By Geroch's theorem, for DFW this is the "spin 
connection" acting on the trivial bundle V x C^. It is built from the "spin 
matrices" F^, [H [21 [3] , 

D^i; = d^^ + T^ij. (5) 
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The formula above may be used also for TRD-1 and for TRD-2, by defin- 
ing appropriate matrices from the relevant connection [22j- The spin 
matrices F^, of DFW, depend on the set of fields (7^^) (/x = 0, ...,3). The 
F^ matrices of TRD-1, and those of TRD-2, do not depend on the fields (7'^). 

The hermitizing matrix is a nonzero 4x4 complex matrix A such that 

[2Sll2n] 

A^ = A, {AYy = Aj^ /i = 0,...,3, (6) 

where = M*-^ denotes the Hermitian conjugate of a matrix M. Previous 
work [TTl [22] has proved the existence and uniqueness, up to a real factor 
A(X) 0, of A = A{X), in any spacetime. It has also proved that, if goo > 
and the 3x3 matrix (gjk) {j,k = 1,2,3) is negative definite, then the sign 
of the matrix A can be chosen such that B = ^7° is a positive definite 
matrix, which is hermitizing for the alpha matrices: 

B^ = B, {Ba^')'^ = Ba" /i = 0,...,3, (7) 

where 

a° = ^Vr, = 7°7V^7°°- (8) 

For DFW, the gamma field "^^{X) is given by Eq. ([1]). It follows easily from 
this ^22j that, for DFW, the hermitizing matrix field can be chosen to be: 

VX G V, A{X) = A\ (9) 

with a constant hermitizing matrix for the "fiat" constant Dirac matrices 
7*^" in Eq. ([T]). Moreover, the standard choices, i.e., the Dirac, chiral, and 
Majorana representations of the constant Dirac matrices 7^" in Eq. ([T]), are 
such that 7''° turns out to be a hermitizing matrix for the set (7^°). Thus, 
in most of the literature, the hermitizing matrix field A is just the constant 
matrix 7^*^, hence one does not usually refer to refer to the notion that the 
constant hermitizing matrix A for DFW can be different from 7^". (Except 
forRefs. PEO].) 

The Dirac Hamiltonian H, which puts the Dirac equation ([3]) into Schrodinger 
form: 

i^ = H^, = (10) 

^ Thus, the hermitizing matrix A{X) is unique up to a positive factor X{X) > 0. 
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is PES]: 

H = m«° - iaWj - iTo. (11) 

[Here Tq = Dq — do is the connection matrix of Eq. ([5]) .] It is important to 
note p2] that, in a given spacetime {y-,g^u) and with given coefficient fields 
(7^^,^), this operator still depends on the coordinate system, or, more ex- 
actly, on the reference frame — the latter being understood here as an equiva- 
lence class F of local coordinate systems (charts) on the spacetime V, modulo 
the purely spatial transformations 

x'O = x\ x'^ = f^iix'')). (12) 

As shown in Ref. [22], the data of a reference frame F determines a three- 
dimensional "space" manifold M, which is the set of the world lines of the 
observers bound to F — i.e., a congruence of time-like curves with spatial 
coordinates x^ parameterizing the space manifold M and a time coordinate 
x° parameterizing the curves. In the present work, we fix the reference frame 
F and even the chart (the latter being, however, an arbitrary one), and 
we study the dependence of the Hamiltonian and energy operators on the 
coefficient fields (7^,74). We give little detail about the covariance of the 
definitions introduced under a change of chart, as this covariance is clear, in 
each case (DFW or TRD), from previous work pTl [22] . 

2.2 Current conservation: the modified Dirac equation 

The probability current is defined as: 

J^' = ip^AY^- (13) 

This definition is generally-covariant, the current being indeed a four- vector, 
for TRD and for DFW as well. The current 0131) is invariant by a change of 
the fields 7^ through a local similarity transformation S = S{X) G GL(4, C). 
This is a transformation that changes one set (7^^) to another one (7''), also 
satisfying (jl]), by setting: 1^1 

7/^ = ^-V^, /i = 0,...,3. (14) 

^ For both the DFW and TRD theories in a non-simply connected spacetime, there 
may be distinct classes of smooth Dirac matrices 7** which do not transform into each 
other by smooth local similarity transformations 26 . However, in this paper, the only 
changes of the Dirac matrices 7^ that we consider are induced from smooth local similarity 
transformations. This restricts the changes of the matrices to be within a single class 
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[Of course, the anticommutation relation (jl]) has to be satisfied by the new 
gamma matrices with the same given metric g^^, as with the starting ones 
j^. It is easy to check this from Eq. f|T^ .] A hermitizing matrix field for the 
new fields 7^ may be defined thus [17J: 

A = S^AS. (15) 
With the corresponding change of the wave function: 

7/- = 5"V, (16) 
this leads to the invariant relation [T71 [22] 

j'' = 4j^ Ay4' = J^- (17) 

Note that the original Dirac theory [TSl [12] employed Lorentz coordinate 
transformations and unitary similarity transformations 5*. The DFW theory 
employs general coordinate transformations and local spin transformations. 
That is to say, local similarity transformations S that, at each point X of 
the spacetime, belong to the spin group Spin(l,3) [a relevant definition of 
it is given after Eq. f l22p below]. The elements of Spin(l,3) are generally 
not unitary matrices. The DFW theory is therefore incomplete for quantum 
mechanics in some aspects by requiring separate treatment for the unitar- 
ily equivalent Dirac representation, chiral representation, and the Majorana 
representation, none of which are related by the spin group Spin(l,3). We 
shall base our discussion of the DFW version of the Dirac equation on the 
following three properties, which characterize the DFW theory as it is used 
in most of the literature: |j 

a) The Dirac matrices 7^^ for DFW are of the form 7^^ = a^^ 7^°", where 
a'^Q, are real coefficients such that a^^ a" ^ i]"^ = g^'^, and where 7^° are a fixed 
choice of constant Dirac matrices, taken from a standard representation {i.e., 
chosen to be in the Dirac representation, or in the chiral representation, or in 
the Majorana representation, or in any other representation that is unitarily 
related to these standard representations.) 

b) The hermitizing matrix for the Dirac matrices 7^^ is A = 7^°. This 
choice of constant hermitizing matrix is invariant under all local spin trans- 
formations S. 



^ See the exceptions mentioned after Eq. 
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c) The coefficient fields (7^, A) for DFW are covariantly constant witli re- 
spect to tlie spin connection of the form: T^, = u^prS^'^ with real coefficients 
Uiypr and sP'^ = [7^'', 7''^]. That is, all covariant derivatives of the coefficient 
fields (7^, A) for DFW vanish. 



In a previous work [22], it was proved that, in order that any solution i/j 
of the Dirac equation ([3]) satisfy the current conservation 

D^J'' = 0, (18) 

it is necessary and sufficient that 

D^B^" = (5^' = AY)- (19) 

Therefore, not all possible coefficient fields (7^^, A) are physically admissible 
for the "normal" Dirac equation but merely the ones which, in addition 
to the anticommutation relation (jl]), satisfy condition f ll9p . Such coefficient 
fields we call admissible for the "normal" Dirac equation ([3]). However, it was 
also proved |22] that the current conservation f lTSj) applies, for any coefficient 
fields, to any ip solution of the following "modified Dirac equation:" 

YDpij = -%mi) - ]^A-\D^B^)^. (20) 

In other words: for the latter equation, all coefficient fields (7^, A), whose 7^^ 
fields obey (jl]), are admissible. Note that the modified Dirac equation ^20^ 
coincides with the normal one ^ for all ip, iff the condition l[W\) is satis- 
fied. In particular, for DFW, the hermitizing matrix field is given by Eq. 
(Q. Then, since the coefficient fields (7^^,74) are covariantly constant, the 
condition f lT9]) is automatically satisfied |22]. Thus, for DFW, the modified 
Dirac equation l[2U\) coincides with the normal one 

However, for TRD theory, the condition for current conservation f[T^ rep- 
resents a nontrivial linear constraint of sixteen scalar PDE's to be satisfied 
by the coefficient fields (7^, A). By considering as the unknown the matrix 
S" of a similarity (!T^ )-( !T5]) . one may go from any possible coefficients fields 
(7^, A) to ones (7^^, A) which are admissible for the "normal" Dirac equa- 
tion on]), by solving a linear system of sixteen scalar PDE's for the sixteen 
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unknowns of S [22]. But, in practice, this would involve heavy numerical 
calculations. Even in the spherical static situation with diagonal, spatially- 
isotropic metric: (g^u) = diag((i^) with do = ip~^, di = d2 = = —i^, 
(f = (p{r), r = ^/ix^)"^ + (x^)^ + (x^)^, we found that the condition ( fT9l) is 
not satisfied by the most obvious choice of Dirac matrices. Namely, for TRD, 
while there exist fields of Dirac matrices 7^ that satisfy the condition (fT9|l . 
they are generally not so readily obtained as in Eq. ([1]). 

Therefore, to use the TRD theory, we must have recourse to the modified 
Dirac equation ( 12 Up with automatic current conservation. At the same time, 
switching from the normal Dirac equation (|3]) to the modified one fl20|) is 
entirely vacuous for DFW, since, as stated above, ([3]) and f l20|) coincide for 
DFW. 

2.3 The Lagrangian and the spin group 

The modified Dirac equation f l20|) is derived from a Lagrangian which gener- 
alizes the well known "Dirac Lagrangian" valid for the DFW equation (e.g. 
[H ED 123] )■ The Lagrangian density is given by: 

/ = ^ [^Y{D,4^) - (^) Yi^ + 2^m^^/;] , (21) 

where ip = ip^A is the generalized adjoint of ip and similarly D^^ip = {Dfj_ip)^ A. 
In the general case that we are considering, it is straightforward to check that 
the Euler-Lagrange equations of this Lagrangian give the generalized Dirac 
equation fl20l) . 

The point dependent gauge transformations that leave the Lagrangian 
invariant form a transformation group, say Ggaugc- This transformation group 
can equivalently be defined as the set of those local similarity transformations 
under which the Dirac equation is covariant. For the DFW theory, the Dirac 
matrices are defined by Eq. ([1]). Then the group Ggauge is made of those local 
similarities which are deduced from a change of the tetrad field Ua = a^a ^m' 
entering Eq. ([T]), by a (proper) local Lorentz transform L = L{X) G S0(1, 3). 
This is well known, see e.g. on p. 468 of Ref. [Ij. Then one may show [22] 
that we have Eq. f|T^ with 

S{X) ^ ±S(L(X)), (22) 
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where L i— )■ ±S(L) is the spinor representation (defined only up to a sign). 
The set of all S(L) for L e S0(l,3) makes the spin group G = Spin(l,3). 
Thus, Ggauge for DFW is made of all local similarities which at any point 
belong to the spin group Spin(l, 3). On the other hand, for the TRD theory, 
the transformation group Ggauge was described by a partial differential equa- 
tion in Eq. (69) of Ref. [22]. 

In Subsect. 13.11 we consider a possibly larger group of local similarity 
transformations which preserve the form of each theory, but are not covariant 
transformations of that theory. 

2.4 Hamiltonian of the modified Dirac equation and 
Fundamental Formula 

Multiplying the modified Dirac equation fl2U]) by on the left, using the 
anticommutation relation (jl]) and the definition ([8]) of the alpha matrices, the 
Hamiltonian operator H' associated with the modified Dirac equation (120|) . 
analogous to the Hamiltonian operator H in Eq. (fTTI) . is given by: 

H' = m«° - laWj -rV^-'-^A-^ {D,,B^). (23) 

2 9 

Using again (jl]), the inverse of the hermitizing matrix B 
to be 



hence we may rewrite 

H' = ma° - laWj - iTq - -B-^D^B^") = H + zC, C = --B-^D^B^"). 

^ ^ (25) 

We will now identify the relevant Hilbert space scalar product between wave 
functions, say {ip \ </?), and derive the hermiticity condition of H' for that 
scalar product. Due to Eq. ( |25l) . these two steps are only a slight modification 
of the corresponding steps |22j for the Hamiltonian H of the "normal" Dirac 
equation, hence we shall be brief. To ensure compatibility with the unique 
scalar product identified for affine coordinates in a fiat spacetime [T7], the 



B = A'y is seen 
(24) 
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scalar product should have the general form[f| 

(V' I (^) = [ {iIj{x) : (^(x)) dV(x) = / tp^B^ a d=^x, (26) 

where {u : f ) is a positive-definite Hermitian product between 4-vectors (for 
TRD) or quadruplets of scalar fields (for DFW), defined by pTl [22] : 

{u:v) = {u: v)b = Bp^u''*v'' = u^Bv, B = ^7°, (27) 

and where a{t,x) is a priori any spatial scalar field [i.e., scalar under trans- 
formations ([T2|) ]. From f l25|) and f l27|) . we find that 

{il) ■ HV) - (HV : ^/;) = (^ : H^/-) - (Ht/' : V) + i[(^ : C^-) + (C^/- : ^)], (28) 

with 

(V; : CiP) = -^ij^BB-\D,B'')ij = -l^tp^^M)^ = ^cij : V^). (29) 
Therefore, using Eq. (92) in Ref. which applies to a time-independent 

{i) : HV^) - (HV^ : V^) = A^i'' D + (Z}^V^)tA7'^V^], (30) 

we get for a such ip: 

{tP : HV) - (HV : V^) = -^ [^'^^t'^^mV' + {D^i))^AYij + 7/>t(Z}^5^)V'] 

= -zZ}^(^U7'^^) 

= -zD^J^ (31) 

This is Eq. (93) of Ref. [22], which was derived there for the normal Dirac 
Hamiltonian f lTT]) . by assuming then the condition for current conservation 
fin?]) . Thus, it is true for the modified Dirac equation, without assuming the 
validity of ffH?]) . By integrating fl^ in space, using the general form fl2Bl) of 
the scalar product, we derive first in exactly the same way as in Ref. [22] : 

I HV) - (HV I ^)] = / V'^ [or 9o (v^ ^7°) - (9,-or) ^ A^^] ^ d^x. 

(32) 

^ The second equality assumes that the "space" manifold M (see the end of Subsect. 
12. ip is difFeomorphic to R'^. This is by no means an essential assumption. 
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From this, follow the same consequences, now for the modified Dirac equa- 
tion fl20l) . as the ones derived in Ref. P2] for the normal Dirac equation 
Namely: 

i) Scalar product. (Theorem 5 of Ref. [22|.) In order that the following two 
axioms be simultaneously valid: 

Axiom (A): The Hilbert space scalar product [2D\] of any time-independent 
wave functions ip and defined on the space manifold M is time independent, 

Axiom (B): For each time t, the Hamiltonian H' is a Hermitian (i.e., 
symmetric ) operator with respect to the scalar product 

it is necessary that the scalar field a be a constant, and thus may be assumed 
to be cr = 1, so that the scalar product must be 

f iiJ: v)b dV = /" ^^B^ d^x, dV = d^x. (33) 
Jm J 

As noted in Ref. [22], Axioms (A) and (B) are a somewhat weaker version 
of two axioms valid for standard quantum mechanics in a flat spacetime. 

ii) Hermiticity condition. (Theorem 6 of Ref. [22j, now without the necessity 
to fulfil the D^B^ = condition.) In order that the modified Hamiltonian 
( 125|) be Hermitian (at time t) for the scalar product (133|1 . it is necessary and 
sufficient that 

aoM(x) = for every x G R^ M = B. (34) 

Moreover, since the scalar product must be (l33i) . i.e., a = 1, Eq. (l32il 
becomes the following Fundamental Formula: 

I HV) - (HV \ij)]= [ tp^ (doM) i/j d=^x. (35) 

This fundamental formula, as well as the condition a = 1 of Eq. ( I32|) . applies 
of course also to the "normal" Dirac Hamiltonian H in Eq. (ITT]) . 
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3 Equivalence condition of the Hamiltonians 
corresponding with two coefficient fields 



3.1 Definition of equivalent operators 

As was noted in the llntroductiont there is an infinity of possible choices for 
the coefficient fields of the Dirac equation in a given spacetime {Y^g^^). As 
discussed in Subsect. 12.21 if (7^, A) and (7^^, A) are two possible choices for 
the coefficient fields, we assume that they are related together by a local 
similarity S{X) as in Eqs. and f|T5l) . (See Footnote 3 there.) We inves- 
tigate the changes in the Dirac Hamiltonian fl25l) . in the energy operator fl60l) 
below, and in its spectrum, that follow from a change of the coefficient fields 
(7^^,74) through a local similarity ffM|) - ffT^ . When doing so, the local simi- 
larity transformations do not necessarily leave the Dirac equation covariant 
under the transformation (fT6l) of the wave function. Indeed, while the Dirac 
equations obtained from different choices of the coefficient fields are not all 
equivalent as wave equations, each of them is derived from a Lagrangian 
and has a conserved probability current. [For the modified Dirac equation 
(1201) . any solution of the Dirac equation with the new coefficients satisfies 
the current conservation (|T8l) . see Subsect. 12.21 .] While the admissible lo- 
cal similarities 5* ensure that the anticommutation relation (jl]) holds true, 
they must also be consistent with the form assumed for the coefficient fields. 
Recall the three properties (a), (b), (c) of the DFW theory discussed in the 
paragraph immediately following Eq. f lT7|) . 

Thus, for the TRD theory, any local similarity S which at any point be- 
longs to GL(4, C) is admissible in this generalized sense [with the modified 
Dirac equation f l20p ]. However, for the DFW theory, one assumes the par- 
ticular form ([T]) for the field 7^^. As recalled in Subsect. 12. 3[ this implies 
that the admissible local similarities belong at any point to the spin group 
Spin(l, 3), thus any of them leaves the Lagrangian fl2T|) invariant. 

With each of the two coefficient fields: [j^jA) and (7'^, A), corresponds 
a unique scalar product: fl55]) for the first one, and 

{x]uj) = J{x:io)^dy = jx^Bu^d?^, B = Af (36) 

for the second one. These two scalar products are isometrically equivalent 
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through the transformation fll6p of the wave function: from fll4p . fllSp and 
(|5^ 9. we get 

B = S^BS, (37) 

whence 

(^1^)^ [ {S-'ijY S^BS iS-'cp)y^d'^ = \ if). (38) 

The Dirac Hamiltonian is H' [Eq. fl25p ] for the coefficient fields (7^, A), and 
its rewriting with tildes for (7^^, A): 

H' = mS-^a°S - iS'^a^SDj - zTq - ^-B-\D^B^). (39) 

For TRD, we have = and = F^. However, for DFW, the covari- 
ant derivatives and the associated connection matrices F^ depend on the 
fields 7^^, hence are denoted by and F^ for (7^^, A). Recall also that, for 
DFW, the last term in Eq. ( 139|1 vanishes, so that the prime may be omitted. 
The general reasoning below applies to any quantum-mechanical operator, 
thus it applies to H and to H' as well; so we shall omit the prime also for TRD. 

The Hamiltonian operator H is acting on a space of functions (or "states" ) 
ip, each of which is defined on the "space" manifold M that is associated with 
the reference frame considered (see the end of Subsect. 12. ip . The set 7/ of 
all wave functions ip defined on M that are square-integrable with respect to 
the scalar product f l5^ form a Hilbert space H. As we know, the operator 
H acts on a dense subspace V of Ti, the domain of H. We shall not need a 
precise definition of this subspace of "H for the operator (125|) . but only the 
fact that it is dense. Being an operator with dense domain "D, H is fully 
determined by the set of the products (Rijj \ ip),ip,ip & V. In particular, 
these products determine the spectrum of the operator H. [The spectrum of 
H is real when H is Hermitian (more precisely self- adjoint). When that is not 
the case, the energy operator E should be substituted for the Hamiltonian 
operator H in all allusions to the spectrum or to the stationary states: see the 
beginning of Section H] for the definition of E and its relation with H.] In the 
same way, of course, the operator H, corresponding just to another possible 



^ Indeed, fix any ip G T), and set 9 = Hip E H. The linear form (p i~-> F{ip) = {0 \ (p) 
is continuous on H. Since V is dense in H, this linear form, hence the vector 9 E %, are 
determined by the values F{ip) for ip E V, thus by the products (H?/) | ip), ip G T>. 
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choice for the coefficient fields, is fully determined by the set of products 
(Hx I ^ ^) with V the dense domain of H in the Hilbert space ?i, 

endowed with the scalar product fl36|) . Therefore, from the isometry fl38|) of 
T-L onto l-L: if) ^ ip^ii follows that H and H are mathematically and physically 
equivalent iff, for all ip^ip in the domain V of H, we have 

(Ht/; I ^) = (HV^ I (40) 

In particular, the operators H and H have then the same spectrum. However, 
we get directly from fl38|l : 

(H^i^) = (HV^|y.). (41) 

Hence, in order that (110|) be true, it is necessary and sufficient that, for all 
eV, we have 

HV = H?/;, (42) 
and since, from flTB]) . we have YLtp = S^^ Hip = S^^ HSip, this rewrites as 

H = 5-iHS. (43) 

This is the condition of equivalence of the Dirac Hamiltonians associated 
with two different choices of the coefficient fields. This notion of equivalence 
is physically relevant to the problems related with the spectrum and the sta- 
tionary states. Note that stationary solutions in the standard sense exist 
for time-independent coefficient fields in a stationary metric. In the general 
time- dependent case, in particular for a non- stationary metric g^^, at each 
given time to, one may "freeze" the Hamiltonian operator Hjg of that time, 
and search for stationary solutions of the corresponding Schrodinger equa- 
tion. One thus obtains, at each time to, an energy spectrum and a set Sf^ 
of "stationary states" . The significance of stationary states, even for a time- 
varying Hamiltonian operator H^^, is that according to quantum mechanics: 
if a particle is in a state ip, then an energy measurement at time t will yield 
an eigenvalue of Ht^. Immediately following the measurement, the state ip 
will collapse into an eigenvector of H^^, as a result of the measurement. 

3.2 Equivalence condition of the Hamiltonians: DFW 

Let us investigate when exactly does a change of the coefficient fields lead 
to an equivalent Hamiltonian, i.e., when do we have Eq. fH51) ? For the 
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standard version of the Dirac equation in a curved spacetime, which we 
denote by DFW, we may consider the "normal" Dirac equation ([3]), hence 
the "normal" Dirac Hamiltonian (11 II) . as noted at the end of Subsect. 12.21 
For DFW, the connection matrices change after the similarity transformation 
([T4|) in the following way [3]: 

t, = S-'T,S + S'\d,S). (44) 

For both DFW and TRD, this is in fact the characteristic relation to be 
satisfied by the 's, in order that the Dirac equation ([3]) be covariant 
under the similarity (fT4|) - (fT6|) — as one easily checks. For DFW, it is always 
satisfied when S has the required form (l22l) . For TRD, the local similarity 
S is arbitrary, but the F^ 's are invariant, so that the Dirac equation 
generally is not covariant under the similarity f lT4|) - f[T6|) [22]. From ( 1TT|) with 
the tilded fields, and from (jS]), we get for DFW: 

H = mS-^a^S-'iS-^a^S{dj + S-^rjS + S-\djS))-tS-^roS-iS-\doS), 
while we get directly from f lTT]) : 

S'-^H^ = mS~^a^S - tS'^a^dj + F^O^ - iS-^ToS. (45) 
Moreover, we have by the composition rule of operators and Leibniz' rule: 
Sdj =djS -{djS), {djS):'ilj^{djS)^lj, (46) 

whence 

H - ^-^H^ = -iS-^a^Sdj - iS~^a^{djS) - iS-\doS) + tS~^a^djS 

= -tS-\doS). (47) 

Therefore, the DFW Hamiltonian operator H corresponding with the new 
(tilded) coefficient fields is physically equivalent to the one, H, corresponding 
with the starting coefficient fields (H = 5'^^HS'), iff 

doS = 0, (48) 

that is, iff the similarity S is independent of the time t = — in the refer- 
ence frame F which is considered, see the end of Subsect. 12.11 Note that, in 
the derivation of the condition for equivalence, Eq. (jUD, we used only the 
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easily checked expression fill I) of the DFW Hamiltonian, together with the 
standard transformation fl44p of the spin connection. 

However, nothing prevents one from changing the coefficients by a time- 
dependent similarity, which leads to inequivalent Hamiltonians. Thus, the 
DFW Hamiltonian operator is non-unique. If the metric g^i, is time- independent 
(in the frame F), we may choose coefficient fields (7^, A) that are indepen- 
dent of t. However, in the general case {g/^ufi 7^ 0), any possible coefficient 
fields will depend on t. In that case, there seems to be no way to define a 
unique Hamiltonian operator. 



3.3 Equivalence condition of the Hamiltonians: TRD 

For TRD, the Hamiltonian with the new coefficients, Eq. ( 139|) . is 

H' = mS-^a'^S - iS-^a^SDj - iTo - ^-B-\D^B>'). (49) 

In TRD, the components S^^^ of the similarity matrix S build a (J) tensor, 
whose covariant derivative may be written in matrix form as 

D^S = d^S + T^S - ST^. (50) 

We have then by the composition rule of operators and Leibniz' rule: 

SDj = DjS - (DjS), (DjS) : ip ^ {DjS)ip, (51) 

and also, since, from and 1^, B^" = S^B^^S and B = B^: 

]3~\D^B^') = S'^B-\S^)-^[S\D^B^')S +{Df,SyB>'S + S^B^'iDf.S)] {52) 
= S-^[B-\D^B^')]S + S-^B-^[B^'{D^S) + {S^y^D^S^ B^'S]. 

So that [using the fact that B^^B^ = 1 and that, from Eq. above, 

B-^B^ = a^: 

H' = mS-^a^S - tS-^a^{DjS - [D^S)) - zFq - ^-S-^[B-\D^,B^')]S 

-'-S-'[iDoS) + a^{DjS) + B-\S^)-\D^SyB^^S]. (53) 
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On the other hand, we have directly from f l25p : 

S-^E'S = mS-^a'^S - iS'^aW^S - iS-^TqS - ^-S-^[B-\D^B>')]S. (54) 

Therefore [multiplying fl53|) and fl54l) by S on the left], H' = S^^H'S* is 
equivalent to 

ia^{DjS) - iSTo - ^[(DoS) + a^{DjS) + B-\S^)-\D^Sy B^'S] = -iro5,(55) 
or [using (IHJ]) ] to 

lla'(D,S) - B-^(Sr WB^S] = l(D„S) + ST„ - T„S . (6„S) - 

which, using the remark made before Eq. (l53ll . may be finally rewritten as 

2B{doS)S-^ = B^{Df,S)S-^ - [B^{Df,S)S-^]l (57) 

The r.h.s. is antihermitian, but not the l.h.s. These are sixteen, generally 
independent, quasilinear PDE's, for the sixteen independent unknowns S^^. 
One may choose freely the field of Dirac matrices 'j^ satisfying the anti- 
commutation (jl]), plus the arbitrary scalar of the corresponding hermitizing 
matrix A{X), thus getting a starting set of coefficient fields, (7'^, A), and 
the corresponding Hamiltonian H', Eq. ([25D- Then, those coefficient fields 
(7^^, A), whose associated Hamiltonian H' is equivalent to H', are likely to be 
determined uniquely, modulo an initial data, by Eq. (l57j) . But the possible 
coefficient fields are deduced from (7^^, A) by a fully arbitrary local similarity 
S{X), which, most generally, will thus not obey Eq. ( 157|) : the vast majority 
of them are hence not equivalent to (7^^, A). In other words, the Hamiltonian 
of the modified TRD equation is not unique, either. 

4 Equivalence condition of the energy opera- 
tors corresponding with two coefficient fields 

The previous section has shown that, for both the standard (DFW) and our 
alternative (TRD) versions of the gravitational Dirac equation, the Hamilto- 
nian is non-unique: one set of possible coefficient fields (7^^, A) being given. 
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the vast majority of the other ones (7^^, A) lead to a Hamiltonian operator 
which does not satisfy the condition of equivalence With this result, the 
foregoing observation [22] that, for DFW, the hermiticity of the Hamiltonian 
is not stable under admissible similarities, is not surprising any more. It also 
indicates that the spectrum of the Dirac Hamiltonian may be nonunique, too. 
However, since the Dirac Hamiltonian H (or H' for the modified Dirac equa- 
tion) is not always Hermitian, the relevant spectrum is rather the spectrum 
of the Hermitian part of H or H': 

ff^i(H + Ht), (58) 

where the Hermitian conjugate operator H^ is w.r.t. the unique scalar prod- 
uct ( l33l) . We get from this definition, by using ( l35i) and ( l33l) : 

(^|H» = ^[(^|(H + H^)^)] 

= i [(H^ I tA) + (^ I H^)] 

= (^ I HV^) + ^ [(H^ I ^) - (^ I H^)] 



(^ I H^) + ^ y ip\doM)ij d^x (M = y/^B) 



ii; I Hz/;) + / iP^M ( ^M-\doM) ) ipd 



3x 



where since M = y/—g B, 

E = ll + —^B-^do{V^B)=}i\ (60) 
2 V 9 

This result applies to the "normal" Dirac Hamiltonian H, Eq. (fTTj) . and to 
the modified one H' f l25p as well. This Hermitian-symmetrized operator E 
or E' coincides, in the particular case envisaged by Leclerc [23], with what 
one (including him) calls the energy operator, which is derived from the field 
Lagrangian. [The case considered by Leclerc is the DFW equation, with 
the tetrad matrix a^^ in Eq. ([1]) satisfying a^j = for j = 1, 2, 3.] As 
with Leclerc, we can show generally that the expected value of the energy 
operator E' equals the classical field energy as derived from the Lagrangian 
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(^T^. Of course, the condition of equivalence of the energy operators E and E 
corresponding with two coefficient fields (7^^, A) and (7^^, A) is the same as for 
the Hamiltonian operators, Eq. f l40|) [with E or E' in the place of H] . When 
it is true, the operators E and E have, in particular, the same spectrum. As 
with H, (j40|) is equivalent to (H3|) . thus here to 

E = S-^ES. (61) 

In order to see when this is true, we will be able to use the work done 
in Sect. [3] for the Hamiltonian, since, from fl57|) and f pU]) . we get after a 
straightforward computation: 

E - S-^ES = H - S-^RS + ^ [S-^B-\S^)-\doSyBS + 5-^(^0^)] . (62) 
This is true for the modified Dirac equation as well as for the normal one. 

4.1 Equivalence condition of the energy operators: DFW 

For DFW, we have from fliTjl and (162]) : 

E - S'^ES = -iS-\doS) + ^ [S-^B-\S^)-\doSyBS + ^"^(9o^)] 

= ^ [S-^B-\S^)-\doSyBS - S-\doS)] , 
SES-^-E = ^B-^ ^^[B{^oS)S-^y - B{doS)S-^y (63) 

Hence, E = S'~^ES' occurs exactly when 5* satisfies 

B{doS)S-' - [B{doS)S-^y = 2 [B{doS)S-'Y = 0, (64) 

which is the vanishing of an antihermitian matrix, involving sixteen inde- 
pendent real quasilinear equations. Any time-independent similarity is a 
solution of (|64l) . For an arbitrary allowed similarity, S{X) = S{L{X)), one 
has in general E 7^ S~^ES, so that the energy operator of the DFW equation 
is non-unique. Moreover, the difference operator SES^^ — E is actually a 
rather general operator, Eq. fl63|) . hence the spectra of the operators E and 
SES*"^ are likely to be generally different. Section |5] proves that, for any 
arbitrary coefficient fields {'y'^,A), one may choose the similarity S, of the 
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form needed for DFW, so that definitely the spectra of the operators E and 
SES~^ are different. But Eq. (HTi) with E in the place of H shows that, for 
any similarity S, the spectrum of E is the same as that of S~^ES, thus sym- 
metrically the spectrum of SES~^ is just the same as that of E. Therefore, 
the energy spectrum of the DFW equation is non-unique. 

4.2 Equivalence condition of the energy operators: TRD 

For the Hamiltonian of the (modified) TRD equation, we get from Eqs. ( 153|) 
and dSH): 

- (^SH'-H'S) = a^{DjS)-B-\S^)-\D^,SyB>'S-{DoS)-2STo + 2ToS, 

(65) 

hence ( !62l) gives us 

^(SE'-E'S) = a^{DjS) - B-\S^r\D^SyB^'S - (DoS) -2STo + 2ToS 
+B-\S^)-\doSyBS + (doS), (66) 

or, using fl50|) and the remark made before Eq. (|53|) . 

B (^SE'S-^ - E') = z [B^'{D^S)S-^ - B{doS)S-^Y . (67) 

Hence, E' = S^^E'S is equivalent to the vanishing of an antihermitian com- 
plex matrix: 

[B''{D^S)S-^ - B{doS)S-^Y = 0. (68) 

Thus, the invariance condition of the energy operator of the modified TRD 
equation imposes only sixteen independent quasilinear real equations ( |68l) . to 
the sixteen independent complex components of the similarity 5*. However, 
not all similarities S{X) (all of which are admissible here) lead to equivalent 
energy operators. For an arbitrary similarity, again the difference operator 
SE'S"^ — E' is a rather general operator, Eq. f l67|) . so that the spectra of the 
operators E' and ^E'^"^ are generally different. Subsection 15.51 shows that 
S may always be chosen so that they are indeed different. And again the 
spectrum of SE'S~^ is the same as that of E'. Hence, the energy spectrum 
of the modified TRD equation (with automatic current conservation) is non- 
unique, too. 
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5 Non-uniqueness proof of the Dirac energy 
spectrum 

In this section we prove the following result: 

Theorem. Let V be any noncompact four- dimensional spacetime that ad- 
mits a spinor structure, whose Lorentz metric g^i, satisfies goo > 0, and whose 
spatial metric hj^ = ~9jk is positive definite. Assume m 7^ in the Dirac 
equation. Then, for DFW as well as for TRD, for any admissible set of co- 
efficient fields (7^, A), there exists an admissible local similarity S such that 
the spectra of the energy operators E and E before and after the similarity 
are different. Consequently, the spectrum of the energy operator for the free 
Dirac equation is not unique in both DFW and TRD. 

The proof given below assumes that the matrix B = Aj'^ is proportional 
to the identity. [For DFW, this is true if one assumes that the matrix {a^^) 
in Eq. ([1]) is such that a^^- = 0: see after Eq. (1771) .] As discussed after Eq. 
fl77|) . the theorem can also be proved without this technical assumption. 

5.1 Perturbation of an eigenstate by a small similarity 

Let us consider a family of local similarities S depending on a small parameter 
e. Thus, 6S being a possible direction for a perturbation of the starting 
coefficients by a local similarity close to the identity, we have 

S{e,X) = I + e{6S){X) + 0{e^). (69) 

The local similarity transformation S{e) transforms the coefficient fields 
(7^^,^) as in Eqs. f|T^ and f|T5l) . and thereby transforms the Dirac equa- 
tion and its associated energy operator E, as well as the Hilbert space scalar 
product in Eq. fl33|) [see Eq. fl36l) ]. Hence the perturbed energy operator, 
denoted as E{e), is generally an operator in a different Hilbert space than the 
Hilbert space in which the unperturbed energy operator E acts. It is there- 
fore convenient to define another operator E{e), having the same spectrum 
as E(£), by bringing back the perturbed energy operator E{e) to act in the 
same Hilbert space as the unperturbed energy operator E as follows: 

E{e) = S{e)E{e)S{e)-\ (70) 
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Let 'ipi^) be an eigenfunction of the operator E(e), with eigenvalue A(e). 
There is a corresponding eigenfunction ipie) = S{e)^^ilj{e) of E(e) with the 
same eigenvalue \{e). Then, using standard perturbation theory of quantum 
mechanics, we may write: 

6X = Xie) - A(0) = (V^(O) I 6E{e)^{0)) + 0(6^), (71) 

with 

SE{£) = E{£) - E(0) = S{e)E{£)S{e)-^ - E. (72) 

Note that this result uses only the hermiticity of the unperturbed operator 
E. 

5.2 General explicit form of the DFW energy operator 

We have for DFW [Eq. 

6E = SES^^ -E = -tB-^ [B{doS)S-^Y , = ^{Q-Q^). (73) 

Recall that, at any spacetime point X, we defined a Hermitian product of 
complex vectors (p71) . by using the positive definite Hermitian matrix B: 

{u : v)b = u^Bv, B = B{X). (74) 

It is evident that the adjoint of a matrix Q with respect to that Hermitian 
product is 

Qt(^) = B-'Q^B, (75) 

i.e., we have 

(gt(^)« : y)^ = (u : Qv)b. (76) 
Using fl75|) . we may rewrite fl73l) as 

SE = -z [{doS)S-'Y'~''^ , Q'^(^) ^ ^-{Q - Qt(^)). (77) 

In DFW theory it is permitted to choose any tetrad field, whose matrix 
(a^„) satisfies Eq. ([2]). It is generally assumed that one tetrad field is not to 
be preferred over another. To demonstrate the non-uniqueness of the DFW 
energy operator it suffices to consider matrices (a'^a) such that a° ■ = 00 Such 

^ The restriction a^^ = applies to the starting tetrad field, which defines the starting 
field of Dirac matrices 7^. We do not impose a similar restriction on the new tetrad field, 
nor equivalently on the similarity S, which defines the new matrices 7''. 



25 



tetrad fields always exist and are almost universally used in the literature. |^ 
Then we get o^q = \^g^ from the orthonormality of the tetrad, hence: 

5 = ^7° = v/?°l4. (78) 

Thus, with the restriction j = 0, which we use here for simplicity, either 
f l73|) or fl77j) is simply: 

5E = -t[idoS)S^']\ (79) 

moreover the product fl7^ is proportional to the canonical Hermitian prod- 
uct in C^, u.v = u^v. However, in the most general case, Eq. (I77|) allows one 
to draw the same conclusions as does fl7^ . 

As is well known: if t S{t) is a smooth curve in whatever matrix group 
G, then ^S~^ belongs to the Lie algebra ^ of G, and conversely any element 
of Q has the form ^S~^ for some smooth curve t ^ S{t) eG.E For DFW, 
S has the form S = ±S(L), where L i— i- ±S(L) is the spinor representation, 
defined (up to a sign) for proper Lorentz transformations: L G S0(l,3), 
Eq. (|22|) . The spin group G = Spin(l,3), set of all such 5* 's, is known 
to be isomorphic to SL(2, C) and hence, as a real manifold, has dimension 
six. The explicit form of the spinor representation [32] is generated by the 
commutators of the "fiat" Dirac matrices 7^": 

s"^ = [7^",7^'^]. (80) 

The six matrices s"'^, < a < /3 < 3, are known to be linearly independent 
(for combinations with real coefficients), thus they form a basis of the real 
vector space Q, the tangent space to G at the identity. Hence, when t h- )■ S{t) 
is any curve in Spin(l, 3), the matrix ^S~^ has the form 

^S-' = u^ps-^ (81) 

^ A version of Cholesky's matrix factorization theorem, which applies to Lorentz metrics 
G = (gfii,) under the assumptions that 500 > and that h = —{gjk) is positive definite, 
states that any matrix (a'^^), satisfying Eq. factors uniquely as a lower triangular 
matrix [t^^), with positive entries on the diagonal, also satisfying Eq. ([2]), and a local 
Lorentz transformation L^^ ^. Thus, in DFW, we can set (a'^^) equal to a lower triangular 
matrix (as in fact it is in most of the literature). 

^ Indeed, ^ is a tangent vector at the point S{t) G G, hence ^S^^ is a tangent vector 
at the identity of G, and thus, by definition, lies in the Lie algebra of G. Clearly, every 
infinitesimal generator of G can be derived this way. 
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for some arbitrary real coefficients cOa/sit) = —upait), determined by the 
arbitrary curve S{t). It follows from this and from (79) that we have 

6E = -i [u^^s'^^Y , (82) 

where, depending on the local similarity S{X) G Spin(l, 3), the six coefficients 
= —uj^a can depend arbitrarily on the spacetime point X. 

From ( IHOl) and the definition of the hermitizing matrix A, we have in 
general 

(s"^)^ = As^'^A-K (83) 
With standard Dirac matrices, for which A = 7^° = A~^, we get thus 

[^a,S-'T= E ^a.is-'-l^'s^-^^'^). (84) 

0<Q!</3<3 

Since the anticommutation relation with the Minkowski metric implies that 

if a = J and /3 = fc, j, k e {1,2,3}, j k, 



otherwise. 
Hence, flM|) is in fact 



(85) 



l<j<fc<3 j,fc=l 



5.3 Specialization to the chiral representation 

For definiteness, we shall now specialize further to the "chiral" Dirac matri- 
ces. However, we mention that, using Eq. f l77|) and some amount of (Dirac) 
algebra, the results can be shown to hold independently of the representation, 
i.e., of the set of "flat" Dirac matrices, and also hold true for any orthonor- 
mal tetrad field a^^. If the 7^" 's are the "chiral" matrices, we get {t32j, Eqs. 
(10.257) and (10.260)}: 

j,fc=i \ / 
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where 9 = {6k) with 6i = 0023 (circular), and where a = (cr^) with a'^ the 
Pauh matrices. Putting together Eqs. ([H2D, dHSD, and (|S7I) . we obtain finally 



6E = P 



N 
N 



) 



= --e.a. 
2 



(88) 



The operator (5E is thus the multiplication by the point- dependent 4x4 ma- 
trix P. Depending on 6, i.e., on the three real numbers ujk, I < j < k < 3, 
the matrix can be in fact any 2x2 Hermitian matrix with zero trace; see 
e.g. Ref. [32], Eq. (5.226). Any such matrix has two opposite real eigen- 
values fi, — /i, and admits an orthogonal basis of eigenvectors, respectively 
M G for /i and v for —fi. Therefore, the 4x4 matrix P has the two same 
opposite real eigenvalues fi, — /x as has A^. 



5.4 Completing the proof for the standard Dirac equa- 



Note that, in the evaluation of 6E (Subsects. 15.21 and l5.3p . we did not use 
the assumption that the local similarity is a small perturbation. Using this 
assumption, we have Eq. (ITTi) . Using the latter, Eq. (ISSil . and the definition 
(155]) of the Hilbert space scalar product, which we rewrite as 



and in which B is given by flTSl) . we get the shift of the eigenvalue as 



where we have decomposed the four-component eigenstate '0 into two two- 
component states: ip = (0;X)- Assume (without any loss of generality) that, 

say. 



tion (DFW) 




(89) 



SX 




(90) 




(91) 
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Select a constant /i > 0. Then, for each x in the space manifold M {x depends 
on the three coordinates x^ , j = 1,2,3), we may select = N{x) in Eq. 
f l88|) such that is an eigenvector of N{x) with eigenvalue /i, so that 

0t]V0 = /i 0t0. (92) 
Since has the eigenvalues fi and —fi, we have then automatically 

X^Nx > -/i xh- (93) 

Here equality is got only if x{^) is an eigenvector of N{x) with eigenvalue 
—fi. In that case (f){x) is orthogonal to x{^) i^i C^. Thus, we conclude that 
> in fl90|) . unless 0(a;) and x(3^) are orthogonal almost everywhere and 
simultaneously / 0V dV = f x^X dV'. 

Already the first feature is impossible if m 7^ 0. It is well known that the 
probability current J'^ satisfies the following Fierz identity [29] : 

J/^J^ = [ip^AY-^] (V^U7^^) = {iP^AiPY + {i-^^ A'j^'^Y. (94) 

The right-hand side vanishes if 0(x) and x{^) are orthogonal in C^. In that 
case, the probability current J'^ lies on the light cone (almost everywhere), 
which is impossible [23] for a solution of a Dirac equation with a nonvanish- 
ing mass term. Q. E. D. 

5.5 The case with TRD equations 

Using (!75)) . the difference operator ( 167)) is rewritten, similarly as we did for 
DFW in Eq. ([77D, as 

6E' = SE'S-^ - E' = iB-^ [B>'{D^S)S-^ - B{doS)S-^Y 

= t[B-'B''{D^S)S-' -{doS)S-'Y'^''\ (95) 

Assuming again for simplicity that the matrix B is proportional to the iden- 
tity, as in Eq. (ITS]) for DFW, this reduces to 

6E' = I {B-^B^[D^,S)S-^ - [dQS)S-^Y = ^ { [i^oS) - (doS) + a^{DjS)] S'^Y . 

(96) 
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The r.h.s. of f l96p shows that the operator SE' is multiphcation by a 4 x 4 
Hermitian matrix which we may denote as P. For TRD, we do not have any 
restriction on the similarity S. We can set 6E' = P = iQ"^ in Eq. (I96l) . Then 
the matrix Q can be chosen freely, by taking the similarity S to solve the 
hnear PDE 

DqS - OqS + aWjS = QS. (97) 

Hence, the Hermitian matrix P = iQ"' can be chosen arbitrarily. Then, the 
reasoning is simpler than for DFW: we get again > by choosing, for 
each X, the matrix P{x) in such a way that the eigenfunction ip oiE takes 
a value %l){x) which is an eigenvector of P{x) for always the same eigenvalue 
/X > 0, for any x G M. In fact, in that case, the unperturbed eigenfunction 
ip is an eigenfunction of the perturbed energy operator, with the eigenvalue 
A + /i. 

6 Conclusion 

In the llntroductiont we asked three questions: i) When does a change of 
the coefficient fields lead to an equivalent Hamiltonian operator? Or ii) to 
an equivalent energy operator? iii) Is the energy spectrum unique? [Note 
that, in experiments, observing different spectra is easier then determining 
whether two energy operators are equivalent, so that (iii) is stronger than 

(ii)-] 

i) In Section [3] we showed, for DFW, that the Hamiltonian operators H and 
H, before and after application of a local similarity 5*, are equivalent iff S 
is independent of time. For the two versions of TRD, the conditions are 
less simple, but only very special local similarities will satisfy them. Thus, 
for both DFW and TRD, the Hamiltonian operator associated with the free 
Dirac equation in a curved spacetime is not unique. 

ii) Similar results were obtained in Section H] for the energy operator which, 
unlike the Hamiltonian operator, is always Hermitian with respect to the 
relevant Hilbert space scalar product as required for an observable. 

iii) For both DFW and TRD, we proved in Section [5] assuming very mild con- 
ditions that, for any admissible set of coefficient fields (7^^,74), there exists 
an admissible local similarity 5* such that the spectra of the energy operators 
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before and after the similarity are different. 

The non-uniqueness of the energy operator is significant because, according to 
quantum mechanics: if a particle is in a state ip, then an energy measurement 
at any given time will yield an eigenvalue of the energy operator (at that time 
in the non- stationary case). Immediately following the measurement, the 
state ip will collapse into an eigenvector of the energy operator at that time, 
as a result of the measurement. Since the Hamiltonian and energy operators 
are not unique, and since there is no preferred way to choose them in either 
DFW or TRD, we cannot at present predict the energy of a free spin-half 
particle in a curved spacetime. It would seem that the only recourse to resolve 
this non-uniqueness of the energy spectrum in a curved spacetime would be to 
experiments. However, the energy shifts due to different "reasonable" choices 
of the coefficient fields (7^^, A) are very small. For example, for rotating Dirac 
matrices in a Minkowski spacetime, a straightforward calculation using Eqs. 
f l82|) and fl86l) shows that the energy shift is of the order of the Mashhoon 
energy term [23]. That means for the current upper limit of angular rotation 
rates of neutron stars, 10^ Hz, the energy shift is two orders of magnitude 
less than the experimental accuracy of the Lamb shift, which is 10^ Hz [19]. 
Nevertheless, in the current state of the theory, nothing prevents one from 
choosing a tetrad field with an arbitrarily large rotation rate, without any 
relation to the rotation of a physical body. 
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